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Abstract 

In this work we study the individual contribution to diffuse 7-ray emission from the secondary 
products in hadronic interactions generated by cosmic rays (CRs), in addition to the contribution 
of vr*^ decay via the decay mode vr*^ 27. For that purpose we employ the Monte Carlo particle 
collision code DPMJET3.04 to determine the multiplicity spectra of various secondary particles 
with 7's as the final decay state, that result from inelastic collisions between cosmic-ray protons 
and Helium nuclei and the interstellar medium with standard composition. By combining the 
simulation results with a parametric model of 7-ray production by cosmic rays with energies below 
a few GeV, where DPMJET appears unreliable, we thus derive an easy-to-use 7-ray production 
matrix for cosmic rays with energies up to about 10 PeV. 

We apply the 7-ray production matrix to the GeV excess in diffuse galactic 7-rays observed by 
EGRET. Although the non-Tr*^ decay components have contributed to the total emission with a 
different spectrum from the vr^-decay component, they are insufficient to explain the GeV excess. 

We also test the hypothesis that the TeV-band 7-ray emission of the shell-type SNR RX J1713- 
3946 observed with HESS is caused by shock-accelerated hadronic cosmic rays. This scenario 
implies a very high efficacy of particle acceleration, so the particle spectrum is expected to con- 
tinuously harden toward high energies on account of cosmic-ray modification of the shock. Using 
the statistics wc find that a continuously softening spectrum is strongly preferred, in contrast 
to expectations. A hardening spectrum has about 1% probability to explain the HESS data, but 
then only if a hard cut-off at 50-100 TeV is imposed on the particle spectrum. 

PACS numbers: 96.50.S-, 96.60.tk, 98.70.Sa, 98.58.Mj 

Keywords: cosmic rays, 7-rays, hadronic interactions, supernova remnants 
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I. INTRODUCTION 



It is generally believed that diffuse galactic 7-rays are a good probe of the production sites 
and also the propagation of accelerated charged particles in the Galactic plane (see 

m 

and references therein). Complementary knowledge of the cosmic-ray (CR) propagation is 
based on the isotopic composition and the content of secondary particles/anti-particles of 
the locally observed CRs which, however, requires a number of model-dependent assump- 
tions. Also, it is not clear whether or not the locally observed CR spectra can be taken as 
representative for the entire Galaxy 

The observed diffuse Galactic 7-ray emission is composed of three spectral components 
0, 0, y, y, 13 • T-rays from bremsstrahlung and from inverse Compton (IC) scattering, and 
also from the decay of secondary particles, mostly neutral pions, produced in collisions of 
CR nuclei with ambient gas. However, the relative contributions of the first two processes 

n 

with respect to the ir^ decays is uncertain [llj], thus leading to a considerable uncertainty in 
a physical interpretation of the observed Galactic 7-ray emission [1^. 

The situation is similarly difficult for individual sources of GeV-to-TeV 7-rays, which 
may have been produced in hadronic interactions or leptonic processes like inverse Compton 



scattering. The most notable examples of sources difficult to interpret are AGN iets 

nn ni — 

1 151 . Ilq and recently observed shell-type SNR's [17|, [181, H 



l4, 



Exact knowledge of the 7-ray source spectrum resulting from hadronic interactions is 
thus of prime importance for a proper physical interpretation. The production of neutral 
pions and their subsequent decay to 7-rays is thought to be the principal mechanism for the 
7-ray hadronic component in cosmic ray interactions, and a number of parameterizations 
have been developed over the years that thus describe 7-ray production in pp-collision 
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20, 121I, I22I, |23|, I2J, |25|. In addition to the neutral pion production, however, the direct 7-ray 



production in cosmic-ray interaction was found important [20], as was the production of 
(S^, S°), {K^, K^), and rj particles, all of which eventually decay into 7-rays [t, 26|. Also, 
the helium nuclei in cosmic rays and the interstellar medium are expected to contribute about 
20-30% to the secondary production in cosmic-ray interactions, for which the multiplicity 
spectra may be different. 

The purpose of this work is to carefully study the 7-ray production from individual decay 
processes, including the direct production, in cosmic-ray interactions. For that purpose we 
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employ the event generator DPMJET-III [27| to simulate secondary productions in both 
p-generated and He-generated interactions. We include all relevant secondary particles with 
7-rays as the final decay products. For the cosmic rays, protons and helium nuclei are taken 
into account. The composition the ISM we assume as 90% protons, 10% helium nuclei, 
0.02% carbon, and 0.04% oxygen. At energies below 20 GeV, where DPMJET appears 
unreliable, we combine the simulation results with a parametric model of vr*^ production 



2J], that includes the resonances A(1232) and A(1600), thus deriving a 7-ray production 



matrix for cosmic rays with energies up to about 10 PeV that can be easily used to interpret 
the spectra of cosmic 7-ray sources. 

As examples we use the 7-ray production matrix to interpret the observed spectra of 
diffuse galactic emission and of one particular supernova remnant. We calculate the spectrum 
of diffuse galactic 7-rays and find, as other groups did before |23|, that the GeV excess is 



probably not the result of an inappropriate model of hadronic 7-ray production. We also 
test the hypothesis that the TeV-band 7-ray emission of the shell-type SNR RX J1713-3946, 
that was observed with HESS, is caused by hadronic cosmic rays that have a spectrum 
according to current theories of cosmic-ray modified shock acceleration. 



II. COSMIC-RAY GENERATED HADRONIC INTERACTIONS 



A. Parameterized vr*^ production in pp collisions 



The 7r° production and subsequent decays are thought to be one of the most principal 7- 
ray production mechanisms in high-energy astrophysics. The vr" production in pp-coUisions 
was well studied both theoretically and experimentally. However, all the current available 



models have met an applicability limit: the isospin model 

■dp 



considering the relation 



'^'^ ^ + —ir- ] 1 is found to be too simplified by comparison with accelerator data 



dp 



28l . l29l . Isd . l3l| for incident proton energies from tens of GeV to a few hundred of GeV, and 



the parameterized vr"*" and tt production cross sections in proton-induced nuclear collisions 
are found to deviate from the experimental data at the incident proton energy Ep > 200 



21 



321]; the scaling model prediction 



GeV or at the transverse momentum px > 2 GeV/c 

[lo| is found to give a poor fit in both differential and total cross sections at low energy 



Ep < 10 GeV [331], and it was also shown 3^ that the kinematic limit of this model is 
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incorrect at energy Ep < 3 GeV; the isobar model [8|, ISJ, l35| is found superior but is valid 
only from threshold up to a few GeV 8j. 

In addition to the theoretical approach, the parametrization of the inclusive cross section 



21 



24|. 



of TT^ production in pp-collisions based on the accelerator data was also studied lol . 
The validity of a parametrization depends on the functional formulae chosen and also the 
data used to determine the free parameters. Thus, the parametrization method often re- 
produces tails extending beyond the energy- momentum conservation limit, which might 
result in unphysical structures when the spectrum is calculated in a broad energy range. It 



is recognized 20|, |24| that these parameterizations lose accuracy in the high energy regime 



and come with a kinematic limit of validity. The parametrization of Stephens and Badhwar 
[lo| under-predicts the yield of high energy pions. The work of Blattnig et al. [21] is valid 
only for energies Ep below 50 GeV, and above this energy it extremely over-predicts the 
pion production and carries a kinematic invalidity; Kamae et al. 2J| have exercised much 
care in determining the inclusive cross section near the pion production threshold, but the 
same problems with kinematic invalidities, etc, arise in the high energy range. 

The spectral distribution of secondary particles and the total cross section are not directly 
obtainable due to the lack of data. Nevertheless, one can calculate the energy spectrum 
or total inclusive cross section by integrating the parameterized differential inclusive cross 
section over the whole phase space (the scattering angle and the energy). By analyzing the 
behavior of the integral, it is possible to evaluate the accuracy and validity limit of different 
parameterizations. For example, by the analysis of the total cross section, a parametrization 
which cannot integrate to a correct total cross section can be ruled out, even though the 
invariant differential cross section data can be reproduced, since the global behavior of 
the parametrization is not accurate. On the other hand, producing a correct total cross 
section upon integration doesn't necessarily imply accuracy of the global behavior of this 
parametrization. This disadvantage results from the fact that the available data are not 
sufficiently constraining the parametrization of differential cross sections. 
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B. Cosmic-ray generated hadronic interactions by Monte Carlo event generator 
DPMJET-III 



In particle physics, the Monte Carlo simulation approach is widely applied to obtain 
information on particle production in hadronic interactions. The principle advantage of 
using the Monte Carlo method is that the codes can be used directly for calculations of all 
secondary products. However, it is quite useful to have simple representations to investigate 
the fundamental spectral characteristics. In this work, the Monte Carlo event generator, 
DPMJET-III, is used to simulate the secondary production in cosmic-ray generated hadronic 
interactions. 

n 

The Monte Carlo event generator DPMJET [27], is constructed for sampling hadron- 
hadron, hadron-nucleus, nucleus-nucleus and neutrino-nucleus interactions from a few GeV 
up to very high energies. DPMJET is an implementation of the two-component Dual Par- 
ton Model for the description of interactions involving nuclei, including photonproduction 
and deep inelastic scattering off nuclei. It is based on the Gribov-Glauber approach and 
treats both soft and hard scattering processes in an unified way. Soft processes are param- 
eterized according to the Regge-phenomenology whereas lowest-order perturbative QCD 
is used to simulate the hard component. Multiple-parton interactions in each individual 
hadron/nucleon/photon-nucleon interaction are described by the PHOJET event generator. 
The fragmentation of parton configurations is treated by the Lund model PYTHIA. 

As a new feature, DPMJET-III allows the simulation of enhanced graph cuts in non- 
diffractive inelastic hadron-nucleus and nucleus-nucleus interactions. For example, in an 
event with two wounded nucleons, the first nucleon might take part in a non-diffr active 
interaction whereas the second one scatters diffractively producing only very few secon- 
daries. Further features of DPMJET-III are a formation- zone intranuclear cascade and the 
implementation of certain baryon stoppings. 

DPMJET-III and earlier versions such as DPMJET-II have been extensively tested 
against data in hadron-nucleus and nucleus-nucleus collisions [271, l36| , for the simulation of 
cosmic-ray air shower problems 37j, and for the atmospheric 7-ray and neutrino production 



381]. Although DPMJET has reproduced well the experimental data in hadronic collisions 



over a quite broad incident energy range, it appears to produce unphysical structure near 
the threshold energy of pion production. Figure [T] shows the vr" and vr^ production cross 
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sections in function of the kinetic energy of protons in pp-co 
III. The experimental values are taken from references fiol . 



lisions simulated by DPMJET- 



39| up to = 2 TeV. To 



be noted from the figure is that the DPMJET results agree with the data except the energy 
range -E^ < 2 ~ 3 GeV. For a high-energy physics event generator this is not uncommon 
behavior. The energy region in question corresponds to the resonance region relating to the 
nuclear cascades in the nuclear interactions. The contribution from the baryon resonances, 
either A(1232) or A(1600), has a very narrow spike-like characteristic at energies close to 



pion production thresho 
differential cross section 



d, and is sometimes added by hand to parameterizations of the 



2^. 



We therefore use the parameterizations 2l|, |2J| including pion production and resonance 



production below a total energy Ep < 20 GeV. For Ep > 20 GeV, the production is given 
by the DPMJET-III simulation; for Ep < 2.5 GeV, it is purely parametric with resonance 
production included; for 2.5 GeV < Ep < 20 GeV, a linear combination in energy is used 
to transition from the parametrization to the DPMJET-III results. The final production 
spectra of secondary particles are smooth and still retain agreement with experimental data. 

We use DPMJET to calculate energy-dependent weight factors, that allow us to paramet- 
rically account for p + ISM and He + ISM collisions in the parametrization approach, that 
has been strictly derived for p + p collisions only. The weights show no strong dependence 
on energy of the projectile particle. 

A remark is in order drawn on the parametrization of pion production, that we use at low 
cosmic-ray energies. The functional form of the pion multiplicity spectra violates kinematic 
limits, so typically a few percent of pions appear to have energies higher than possible. A 
truncation of the multiplicity spectra at the kinematic limit would modify the mean energy 
of pions, thus changing the 7-ray spectra. As the parametrization is a best fit to data, any 
modification would in principle deteriorate the fit. We therefore decided not to correct for 
the violation of kinematic limits. 



C. Decay channels of secondary particles to 7-rays 

The isobar model is based on the theory that tt production in p + p collisions near the 
threshold is mediated by the excitation of the A3/2 isobar, which subsequently decays into 
a nucleon and pions. In this work, the produced isobar is assumed to carry momentum 
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directly forward and backward in the center-of-mass system of the coUisions with an equal 
probability, and then decays isotropically into a nucleon and tt^'s. 

The decay modes and decay fractions for resonances are considered as follows based on 
the Particle Data Group. 



A(1232) - 




A(lDUUj - 


1 


A(1600) - 


A(1232) +7r°; 


A(1600) - 


N(1440) + 7r°; 


N(1440) - 




N(1440) - 


A(1232) + 7r°; 


N(1440) - 


p + 27r°. 



The TT^ contribution from resonances is then added to the non-resonant 7r° production. 
While running the DPMJET-III, all secondary particles produced in simulated hadronic 

interactions are recorded. The following decay modes shows all the decay processes of 
secondary products with 7-rays as the final decay particles that are taken into account in 
the current work. 

• baryonic decays: 

A ^ n + 7r°; 
A ^ n + 7r°; 
E"^ A + 7; 
E+ ^ p + TfO. 

• mesonic decays: 

n° ^ 27; 

7r° — >• e~ + e"*" + 7; 
K+ ^7r+ + 7r°; 
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K- ^ 71- + 7r°; 
Kl ^ 27; 
Kl SvrO; 
Kl^n+ + TT- + 7r°. 

For the two-body decay process, it is easy to evaluate by particle kinematics while for a 
three-body decay process, the Dalitz-plot is used. For the fraction of each individual decay 
process, the values published in Particle Data Group are applied. 

Figure [2] shows the energy spectra of 7-rays contributed from each type of secondary 
particles produced in p+ISM interactions for two different incident proton energies: 150 GeV 
(upper curve) and 10 TeV (lower curve). The ISM composition is taken as 90% protons, 10% 
helium nuclei, 0.02% carbon, and 0.04% oxygens. For ease of comparison, baryonic (left) and 
mesonic (right) decay contributions are separately shown. The total energy spectra of 7-rays 
(the sum of each contribution), the directly produced 7-ray photons in the hadronic collisions 
and the 7-rays from decays are also plotted in each figure for comparison. To be noted in 
the figures is that, in addition to 7r° decays, the directly produced 7-ray photons appear more 
important than other mesons and baryons over almost the whole energy domain. However, 
the secondary mesons contribute more at the peak of the 7-ray spectrum. Integrated over 
the 7-ray spectrum, the other secondary particles have contributed about 20% 7-ray photons 
relative to the vr" decays. 



Some authors 



25 



26l | suggested that a few percent of 7-rays should be contributed by 
decays of rj mesons, but no r] meson is listed in DPMJET read-out table. This is due 
to the fact that with our simulation setup, the r] mesons have already decayed and only 
their decay products are listed. We have verified that the rj mesons are properly accounted 
for by performing test runs under different PYTHIA parameters, so the //-decay is avoided. 
Another independent verification arises from the very small contribution to e^/e^ and fi^ / fi~^ 
pairs from 77 decays, which are found consistent within the statistical uncertainty with all 
77-mesons having decayed. 
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III. SPECTRUM OF 7-RAYS GENERATED BY COSMIC RAYS 



A. Energy loss in production collisions 

The energy loss resulting from the inelastic production collisions can be determined by 
integrating over the yield in all secondary particles. Figure [3] shows the energy loss rate as 
function of the incident particle energy in the p + ISM collisions. To be noted from the 
figure is that the energy loss rate gradually increases with energy above 10 GeV on account 
of an increasing importance of secondary particles other than pions. The thick dashed curve 
in the figure is a functional approximation to the energy loss rate in p + ISM production 
collisions over the entire energy range described as 

- ~ 6.78 X lO'^^^OLZ^ - 1.55 X 10"^^ log(7) + 1.80 x 10"^^ log^(7). (1) 
E 7(7+1) 

To be noted is that the first term in Eq. ([T]) is a functional approximation introduced earlier 
for pion production up to 100 GeV [13]. This term is a good approximation below 100 GeV 
where pions are the dominant secondary product. 



B. The 7-ray production matrix 

For the cosmic-ray hadronic interactions and the subsequent decays, the spectrum of a 
final secondary particle is described by 

«-''^) = I^H^ = 1. '"-(^-«) KlE ) 

with NcR{Ecfi) = ^^^^^y (GeV cm^)~^ defined as the differential density of CR particles 
{p or a). The differential cross section of a secondary particle produced in the {p, a) + ISM 
collision is 

da , , dn . . 

— {Ecn,E) = a,ro,^ (3) 

where (Jprod is the inelastic production cross section, whose value is simulated by DPMJET, 
and ^ is the multiplicity spectrum of the secondary particle in question. 

We can follow the decay of unstable secondary particles into 7-rays, which is a purely 
kinematical problem, and thus obtain the 7-ray spectrum from secondary-particle decay 
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generated by a cosmic-ray particle with energy Ecr as 

Q^(E^) = J2 ^isM [ dEcR Ncr{Ecr) c(3cr cj{Ecr) ^(^7, Ecr) (4) 

where is the 7-ray spectrum resulting from the initial multiplicity spectrum of the 
secondary particle species k or the directly produced 7-ray spectrum in cosmic-ray generated 
interactions. 

If the total energy of cosmic-ray particles {p or a) is parameterized as 

Et = 1.24 ■ (1 + 0.05)^' GeV/n (5) 

and the 7-ray energy is sampled as 

E^ = 0.01 ■ (1.121376)'-°-^ GeV, (6) 

both with sufficient energy resolution, then the production integral, Eq. (j4j), can be written 
as a sum. 

Q,{E,) = Y^njsM^E, NcRiE,)c/3,a{E,) ^ ^(E„E,) (7) 
j k 

= ^nisM^EjNcR{E,)c(3ja,Mij. (8) 
3 

The problem is thus reduced to a matrix operation, in which a vector, that describes the 
cosmic-ray flux at various energies according to Eq. (JSj), is transformed into a vector com- 
posed of the 7-ray source function at a number of energies predefined in Eq. ([H]). 

Eq. ([H]) is the final formula with which one can evaluate the 7-ray spectrum generated by 
cosmic rays. It includes cTj, the production cross section at the cosmic-ray (p or a) energy Ej 
defined in Eq. ([5]), and the 7-ray energy spectrum matrix My. Each element of the matrix, 
Mjj, is the value of the resultant particle energy spectrum j^\E^^Ei,EcR=Ej) with j being the 
index for the generating cosmic-ray particle {p or a) and i being the index indicating the 
7-ray energy. 

In addition to the 7-ray production matrix, we have also produced the production ma- 
trices of all particles that are stable on a time scale relevant for cosmic-ray propagation, 
namely p, p, e=^, z/e, z/e, u^, and v^. 
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C. Examples of 7-ray spectra 



With this production matrix, we calculate the 7-ray spectra for two examples of cosmic- 
ray spectra: a simple power-law spectrum and a power-law spectrum with an energy cutoff: 

NcR cc ^ (9) 
NcR oc Eel''' ■ exp {^^^ , with E,,,,^ = 100 TeV (10) 

V -f^cutoff / 

EcR being the cosmic-ray total energy and both assuming the cosmic-ray composed of 
protons and helium nuclei. 

Figure H] shows the energy spectra of 7-rays contributed by the secondaries produced in 
CR {p and a) + ISM interactions with a simple power-law cosmic-ray spectral distribution 
described in (Q. Figure [5] shows the same as Figure H] but with an energy cutoff in cosmic-ray 
spectrum at E = 100 TeV. The contribution from decays of secondary baryons and mesons 
produced in cosmic-ray generated interactions on ISM are separately plotted for comparison. 
In each figure, the total 7-ray spectrum and the components from 7r° decays and the direct 7 
production are also shown. As shown in the figures, the 7r° decay is the dominant component 
in the diffuse 7-rays spectrum, however, at least about 20% of the 7-ray photons result from 
secondary products other than pions, including the direct 7-ray production. With the energy 
cutoff in cosmic-ray spectrum, it is found that the 7-ray spectrum is modified also at energies 
below 1 TeV, where the spectra in Figure [5] are notably softer than those shown in Fig HI 
This is because the multiplicity spectrum of a secondary particle has a long tail, and a 
high energy cosmic-ray particle can generate several generations of production collisions, 
multiplying the production at lower energies. 



D. The 7-ray spectrum in the GeV band 

The observed GeV-bump in diffuse 7-ray spectrum [gI is a long-term standing mystery in 
7-ray astrophysics. Both hadronic and leptonic scenarios have been proposed to explain the 
GeV excess. A harder spectrum of cosmic-ray nucleons is required for 7-rays of nucleonic 
origin to account for the observed intensity in the multi-GeV range [3] , without violation of 
the observed intensity at a few hundred MeV. However, significant variations of the cosmic- 
ray spectrum in the Galaxy are unlikely jsj . Strong fluctuations in the spectrum of cosmic-ray 
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electrons must be expected, though 0,|4^: Electrons are accelerated by SNR with power-law 
spectra of index -2.0 or similar, but severe energy losses prevent them from homogeneously 
filling the Galaxy. A hard spectral component from inverse Compton scatterings could then 
explain the GeV bump in the diffuse 7-rays. It has also been argued that both a harder 
electron and a harder nucleon spectrum may be required to reproduce the GeV excess jl2 |. 

Figure [6] shows the observed GeV-band 7-ray emission from the inner Galaxy in com- 
parison with the contributions from vr" decay as well as bremsstrahlung emission, which we 
here describe by a power-law spectrum $b(-S) such that 

$b(^) - 1.3 X 10-« ^ . . ( ^ ) """'^ (11) 

^ ^ 0.1 eV/cm^ 1022 cm-2 VlOOMeVy cm2 sec sr ^ ^ 

with power-law spectral index = 2.1, the electron energy density Ue = 

0.1, 0.4, 0.8 eV/cm^ and the gas column density Njsm = 3 x IO22, 8 x 10^^, 3 x 102icm-2, 
respectively. The generating cosmic-rays are assumed with an energy density pe = 
0.75 eV/cm^. Models based on the locally observed cosmic-ray spectra generally predict 
a softer spectrum for the leptonic components, even after accounting for inverse Comp- 
ton emission so we may in fact overestimate the GeV-band intensity of the leptonic 
contribution. Nevertheless, it is clearly seen in this figure, that in the total intensity an 
over-shooting around ~ 300-600 MeV appears in the modelled 7-ray energy distribu- 
tion, whereas a deficit is present above 1 GeV. The observed spectrum of diffuse emission 
is always harder than the model spectrum, and we therefore conclude that an inaccurate 
description of hadronic 7-rays is ruled out as the origin of the GeV excess, in line with earlier 
work 9|. 

IV. THE TEV-BAND SPECTRUM OF RX J1713-3946 

As SNR's are believed to be the main sources of Galactic cosmic rays, it is worth testing 
whether or not the observed 7-ray spectra from SNR can be well modelled in terms of 
hadronic interactions. The 7-ray production matrix presented in this paper is a valuable 
tool in such an analysis. 

The HESS collaboration has measured the TeV-band gamma-ray spectrum of the shell- 
type SNR RX J1713. 7-3946 over more than two decades in energy and in approximately ten 



angular resolution elements per diameter [l9| . The data we use here consist of 25 independent 
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spectral flux measurement for the entire remnant, taken with the complete HESS array in 



2004. Since SNR RX 
along the line of sight 



1713.7-3946 appears closely associated with dense molecular clouds 
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42l | , it is reasonable to expect a strong signal of hadronic 7-rays. 
Given the integrated flux in the TeV band, for the hadronic scenario the required total 
energy in accelerated cosmic-ray nucleons is 

tyr-(2.10»erg)(i^)^(j^^^)"' (12) 

n 

There is kinematic evidence that the distance D is 1.3 ± 0.4 kpc [41], whereas the weakness 
of thermal X-ray emission implies a density of n < 0.05 atoms/cm~^ in the center of the 



remnant [43| and a postshock density of n < 0.02 atoms/cm~^ [4j]. The low postshock 
density indicates that about 1 solar mass or less have been swept up by the SNR, which 
therefore should just be in transition from the adiabatic phase to the Taylor-Sedov phase, 
when its ability to accelerate particles is presumable at its peak. The low density also implies 
a total energy in cosmic rays of about 10^^ erg, somewhat larger than the canonical value 
for the total kinetic energy in a supernova explosion. 

A high energy density of energetic particles is in line with the notion of Fermi-type accel- 
eration at the forward shock of SNR. This process is intrinsically efficient, and a significant 
fraction of the pressure and energy density will be carried by energetic particles. In that 
case the shocks should be strongly modified, because the energetic particles have a smaller 
adiabatic index and a much larger mean free path for scattering than does the quasi-thermal 
plasma. In addition the particles at the highest energy escape, thus making the shock es- 



sentially radiative and increasing the compression ratio 45|. 



While the details of the published acceleration models differ, they essentially all predict a 
continuous hardening of the particle spectrum with increasing energy, until a cut-off occurs 
at an energy that is difficult to predict [46^. At particle energies in excess of 1 TeV, i.e. 



relevant for 7-ray production above 200 GeV, the spectrum should 



display a positive curvature, unless the cut-off has already set 



in H. 



3e harder than E ^ and 



We therefore chose to 



parametrize the spectrum of accelerated hadrons in RX J1713. 7-3946 as 



N{E) = No e [E^ax - E] (13) 

where B is the step function and Eq = 15 TeV is a normalization chosen to render variations 
in the power-law index s statistically independent from the choice of spectral curvature, a. 
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The cut-off energy, i?max, is a free parameter, but should be > 1 PeV, if SNR are to produce 
Galactic cosmic rays up to the knee at a few PeV, because RX J 1713. 7-3946 appears to 
already turn into the Sedov phase. For cosmic-ray modified shocks we would expect s < 2 
and (T > 0. The normalization, Nq, is not a parameter in our analysis, because we normalize 
both the data and the model to their value at 0.97 TeV, while propagating the flux error at 
that energy to the errors in the normalized flux at all other energies. 

The parametrized particle spectrum Eq. f|T3l) is used to calculate the 7-ray source function 
according to Eq. ([8]), which then can be compared to the HESS data. We use the statistic 
to both determine the goodness-of-fit and the confidence ranges of the three parameters, 
-^max, s, and a. The best fit has a good quality with Xmin — 22.48 for 22 degrees of freedom 
and is achieved for s = 2.13, a = —0.25, and -Emax ^ 200 TeV, i.e. no cut-off. The 
best fit is thus not commensurate with expectations based on acceleration at a cosmic-ray 
modified shock. Figure [7] shows the data in comparison with the best-fit model of hadronic 
7-ray production. To be noted from the figure is the obvious value of data in the energy 
range between 1 GeV and 200 GeV as may be provided by GLAST in the near future. 
The presently available upper limits from EGRET 48|] are too high to provide meaningful 
constraints and are therefore omitted here. 

More important than the best fit are the confidence ranges of parameters, for which the 
probability distributions follow a Xn distribution in — Xmin' where n is the number of 
parameters, for which the combined confidence volume is determined j^sj]. As only two- 
dimensional contours can be unambiguously displayed, we have performed a 2-parameter 
analysis by marginalizing one parameter. The areas for the three confidence levels 68%, 
95%, and 99.7%, corresponding to 1, 2, and 3 sigma, respectively, are shown in Figure [8] for 
the curvature a versus the power index s, and in Figure [9] for cr and s versus the cut-off 
energy E^ax- 

To be noted from the figures is that a positive spectral curvature, cr > 0, is excluded at 
more than 95% confidence. This is clearly at odds with expectations based on cosmic-ray 
modified shocks. The contour for -Emax are open toward higher energies, indicating that 
a cut-off is not statistically required. A very high -Emax would be in line with cosmic-ray 
acceleration to the knee at a few PeV, but one should note from Figure [9] that in that case 
a negative spectral curvature is strongly required, in fact a < —0.1 with more than 99.7% 
confidence! A small or marginally positive curvature is permitted only in combination with 
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a low cut-off energy -Emax 100 TeV. 

Tfiese statements are verified by additional tests in which we allowed only a positive 
spectral curvature, u > 0. The best fit under those constraints achieved = 29.39, higher 
than the unconstrained best fit by Ax^ = 6.91, which again demonstrates how much poorer 
the fit is for a > 0. The two parameters -Emax and s would be slightly correlated with pivot 
point at s ~ 1.8 and £^max — 75 TeV, so a lower cut-off energy would compensate for a harder 
power index and vice versa. Nevertheless the cut-off energy would be -Bmax < 140 TeV with 
99.7% confidence. 

These findings clearly demonstrate that a softening particle spectrum is statistically re- 
quired, if the TcV-band 7-ray spectrum from RX J1713. 7-3946 is of hadronic origin. A more 
gradual decay of the spectrum is preferred over a hard cut-off. With the best-fit curvature 
a — —0.25 the particle spectral index would change by As = 0.58 for each decade in energy, 
or by As = 1.04 between our normalization energy Eq — 15 TeV and 1 PeV. There is no 
evidence for either efficient nucleon acceleration to a energies near the location of the knee, 
or for nonlinear acceleration at cosmic-ray modified shocks. 

V. CONCLUSION 

The diffuse galactic 7-ray emission provides information to obtain a proper understanding 
of the origin of galactic cosmic rays. In this work, the 7-ray spectrum generated in cosmic- 
ray interactions with ISM gas is calculated using the simulation code DPMJET-III. For 
the first time we consider both protons and helium nuclei in cosmic rays and include all 
secondary products with 7-rays as the final decay state as well as direct 7-ray production 
and the intermediate resonance production. We find that about 20% of the 7-ray spectrum 
come from non-7r° secondaries, to which directly produced 7-ray photons contribute stronger 
than other secondaries. We also calculate the energy loss rate of hadronic cosmic rays in the 
Galaxy and present a simple analytical approximation to it, that may be used in cosmic-ray 
propagation studies. 

The main product of this work are matrices, or look-up tables, that can be used to 

calculate the source spectra of 7-rays and other secondary products arising from inelastic 
collisions of cosmic rays with arbitrary spectrum. Our study involves 7-rays with energies 
up to 10^ GeV and incoming cosmic rays up to 10* GeV. The matrices will be made available 
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for download at ]http:/ / cherenkov. physics. iastate. edu/ gamma-proa 



We apply the 7-ray production matrices to two examples: the GeV excess and the TeV- 
band 7-ray spectrum of the shell-type SNR RX J1713-3946. Our findings are: 

• The modifications in the GeV-band 7-ray emission of hadronic origin are insufficient 
to explain the GeV excess in diffuse galactic 7-rays. It appears that, to arrive at an 
understanding of the GeV excess, GLAST data are needed to better understand its 
angular distribution and exact spectral form. 

• If the TeV-band spectrum of RX J1713-3946 as observed with HESS is caused by 
cosmic-ray nucleons, then a soft cut-off at about 100 TeV is statistically required in 
the particle spectrum. There is no evidence for efficient nucleon acceleration to energies 
near the knee in the cosmic-ray spectrum. We also find no evidence of the spectral 
curvature and hardness predicted by standard models of cosmic-ray modified shock 
acceleration, even though a strong cosmic-ray modification must be expected given 
the high cosmic-ray density in the remnant implied by the low thermal gas density 
and high 7-ray luminosity. We again emphasize the need for GLAST data to better 
constrain the 7-ray spectrum below 100 GeV. 
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FIG. 1: vr*^ and tt^ production cross sections in pp collisions simulated by DPMJET-III, here shown 
as function of the kinetic energy Ek of incident protons. The data shown figure are quoted from 
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FIG. 2: The energy spectra of the 7-ray contributions from each secondary particle produced in 
p + ISM interactions for two different incident proton energies: 150 GeV (upper panel) and 10 
TeV (lower panel). For ease of comparison, baryonic (left) and mesonic (right) decay contributions 
are separated. In each figure, the energy spectrum of total 7-rays (the sum of all contributions) 
and also the spectrum of the directly produced 7-ray photons are shown. 
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FIG. 3: The energy loss rate for p + ISM collisions as calculated by integrating the secondary 
particle yield. The thick dash-dot line shows an analytical approximation. The ISM number density 
is assumed riisM = 1-0 /cm^. 



22 




FIG. 4: The energy spectra of 7-rays contributed by the secondaries produced in CR (p and a) + 
ISM interactions with a simple power-law cosmic-ray spectral distribution in Figures are 
shown for baryonic (left) and mesonic (right) contributions. The total 7-ray spectrum, the com- 
ponents from TT^ decays and the direct 7 production are also shown in each figure. Note that we 
only sample cosmic rays up to 10^ GeV. 
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FIG. 5: The energy spectra of 7-rays contributed by secondaries produced in CR {p and a) +ISM 
interactions, assuming an exponential cutoff in the cosmic-ray spectrum at 100 TeV described in 
Eq. (jlOp . The plots are arranged as in Figure |H 



24 




llog(E/GeV) 



FIG. 6: Diffuse 7-ray spectrum at GeV range, shown in comparison with data for the inner part of 
Galaxy at 315° < I < 345° and \b\ < 5^ \ The model spectrum is composed of a hadronic com- 
ponent calculated with cosmic-ray spectrum according to Eq. (jlO|) . plus a power-law component 
with spectral index = 2.1 for bremsstrahlung. Different curves correspond to different normal- 
izations of the hadronic and leptonic contributions, here by varying the electron energy density 
uje = 0.1, 0.4, 0.8 eV/cm^ and the gas column density Nism = 3 x 10^^ 8 x 10^^, 3 x lO^^cm"^^ 
The energy density of hadronic cosmic rays is kept constant at pE = 0.75 eV/cm^. 
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FIG. 7: The TeV-band 7-ray spectrum observed from RX J1713-3946 with HESS, shown in com- 
parison with the best-fit model of hadronic gamma-ray production. Note that both data and model 
arc normalized at 0.97 TeV. Data in the energy range between 1 GeV and 200 GeV would be very 
valuable and may be provided by GLAST in the near future. 
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FIG. 8: The confidence regions for the spectral curvature a and the spectral index s. The dark 
shaded area corresponds to 68% probability, or 1 sigma, the medium gray area indicates 95% 
confidence, or 2 sigma, and light gray is for 99.7% probability or 3 sigma. 
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FIG. 9: Upper: The confidence regions for the spectral curvature a and the spectral index s, versus 
the cut-off^ energy i^max- The contour levels are 1 sigma, 2 sigma, and 3 sigma as in Figure [HI 
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